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Abstract — In this research , a study was done on wave profile and type at breaker point. The study was done using 
shoaling-breaking model. Output of the shoaling-breaking model becomes the input at water wave surface equation 
to obtain wave profile. By obtaining the wave profile, the calculation of breaker height becomes more accurate. 
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I. INTRODUCTION 

Wilson (1963) obtained that there are four types of water 
wave, i.e. Airy’s waves, Stoke’s waves, Cnoidal waves and 
Solitary waves. Airy’s type, where wave crest is 
symmetrical with wave trough and wave height is twice of 
that wave amplitude. This type of wave is found only in a 
wave with a very small amplitude. Next is Stoke’s type, 
formed in a wave with wave amplitude bigger than wave 
amplitude of the Airy’s type. In this Stoke's type, the wave 
profile is somewhat asymmetricand the wave height is still 
close to twice the wave amplitude. In a bigger wave 
amplitude, cnoidal type of wave is formed. In cnoidal type, 
the wave profile is asymmetric where the distance between 
wave trough and neutral line or still water level is very 
much different with the distance between wave crest and 
still water level. For a wave with bigger wave amplitude a 
perfect Cnoidal profile will be formed with wave trough 
that almost coincides with still water level. This type of 
wave is also called Solitary wave. Either Cnoidal type or 
Solitary wave type shows that almost all part of the wave is 
above the still water level, nevertheless wave height that is 
twice wave amplitude might occur. 

In shoaling-breaking model, some use wave amplitude as 
the variable and some other use wave height as the variable. 
For this second type, the formulation was done with an 
assumption that wave height has a value twice that of wave 
amplitude. With those four types of wave, then the output of 
shoaling breaking model that uses wave amplitude as the 
variable or the model that uses half of the wave heightcan 
not confirm the wave height. To obtain a certainty of the 
wave height as the result of shoaling-breaking analysis, the 
analysis of water wave surface profile needs to be done to 
obtain a better wave height value. 


This research used water wave surface equation of 
Hutahaean (2019a) to conduct wave profile analysis at 
breaker point. The equation used wave constant G, wave 
number fcand wave amplitude rias the parameter. The three 
wave parameters were obtained from shoaling-breaking 
analysis of Hutahaean (2019b). 

At the same time, this research is also a revision on 
Hutahaean (2019b), where at the research the wave 
amplitude as the result of shoaling-breaking model was 
considered has a value of half wave height. Even though it 
provides a result that is quite close to breaker height from 
breaker height index equation as the result of laboratory 
analysis, the result is less than accurate. 

II. SEVERAL TYPES OF WAVE 

The general form of water wave profile is shown on Fig.l. 
The fonn is assymetrical between wave crest and wave 
trough. If the wave crest elevation against still water level is 
called r ) max ; whereas the wave trough elevation is called 
qmin’ then there is an asymmetry where q max > \q m in I, 
whereas wave height is H = r) max — q min . 



Fig.l. General profile of water wave surface 
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2.1. Wilson Criteria (1963) 

Wilson (1963) classified wave based on the value of ' lm “ x 
(Table (1)). Type Airy’s waves have the value of Vmax < 

H 

0.505, which shows that the wave form is asymmetrical, 
where? 7 ma;l . =» |? 7 min |. The second type is Stoke’s wave 
type that is still quite close with the Airy’s waves 
characters. Furthermore, cnoidal waves type is absolutely 
asymmetrical and the last is solitary waves type, 
where^ 2212 ^ = 1.0. 

H 


Table.1: Wilson Criteria (1963) 


Wave Type 

Vmax 


H 

Airy waves 

< 0.505 

Stoke’s waves 

< 635 

Cnoidal waves 

^ r- Vmax 

0.635 <-< 1 


H 

Solitary waves 

= 1 


In the next sections the analysis of wave profile in the deep 
water is presented, includes the calculation of r ) max , ri min , 

wave height H and wave height ratio jand the depiction of 
the wave profile using equation (1). The result of the 
calculations shows a compatibility with Wilson criteria 
(1963). 

2.2. Airy’s Waves Type 

Airy’s waves type have symmetrical profile where? 7 ma;l . = 

I rimin land wave height H = 2 A. This type of wave is found 
only in a wave with a very small amplitude. Table (2) 
shows the measurement of Airy’s wave profile in the deep 
water for several wave periods with a very small wave 
amplitude A, where = 0.503 was obtained, which 

shows that it is not really symmetrical, but H = 2 A was 
obtained. The profile of Airy’s wave can be seen in Fig. 2. 
Table. 2: Wave amplitude for sinusoidal wave, in the deep 


water 


T 

A 

Vmax 

H 

H 

(sec.) 

(m) 

H 

(m) 

A 

8 

0,017 

0,503 

0,035 

2 

9 

0,02 

0,503 

0,041 

2 

10 

0,023 

0,503 

0,047 

2 

11 

0,027 

0,503 

0,053 

2 

12 

0,03 

0,503 

0,06 

2 

13 

0,033 

0,502 

0,067 

2 

14 

0,037 

0,502 

0,073 

2 


15 

0,04 

0,502 

0,081 

2 



Fig.2. Aiiy’s wave T = 8 sec., A = 0.017 in, in the deep 


water. 


2.2. Type Stoke’s wave. 

Table (3) shows wave profile in the deep water with a quite 
large wave amplitude, with ' ) "^ a * = 0.635that according to 
Wilson criteria is the maximum limit of Stoke’s wave type. 
The value ofj = 2.099, is still quite close to 2. Therefore, 

the characteristics of Stoke’s wave is that wave profile is 
asymmetric but H ~ 2A, although it is at its maximum 
limit. 

Table.3: Wave profile of Stoke’s wave at its maximum limit 


T 

A 

V max 

H 

H 

(sec.) 

(m) 

H 

(m) 

A 

**max 

7 

0,47 

0,635 

0,986 

2,099 

8 

0,614 

0,635 

1,288 

2,099 

9 

0,777 

0,635 

1,631 

2,099 

10 

0,959 

0,635 

2,013 

2,099 

11 

1,161 

0,635 

2,436 

2,099 

12 

1,381 

0,635 

2,899 

2,099 

13 

1,621 

0,635 

3,402 

2,099 

14 

1,88 

0,635 

3,946 

2,099 

15 

2,158 

0,635 

4,529 

2,099 


2.3. Cnoidal Wave Type. 

Table (4) is wave profile in deep water with large amplitude 
where ' lmax = 0.803 was obtained with the value of - = 

H A 

2.726 > 2. Therefore, in this cnoidal wave type the wave 
profile is asymmetric and H =£ 2A. If it is approximated 
with Airy’s waves type there will be a quite large error, in 
this case is 36.2 %. For a wave with bigger wave amplitude, 
there will be even bigger error 
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Table.4: Cnoidal wave profile 


T 

A 

V max 

H 

H 

(sec.) 

(m) 

H 

(m) 

A 

7 

0,899 

0,803 

2,449 

2,726 

8 

1,174 

0,803 

3,199 

2,726 

9 

1,485 

0,803 

4,049 

2,726 

10 

1,834 

0,803 

4,998 

2,726 

11 

2,219 

0,803 

6,048 

2,726 

12 

2,641 

0,803 

7,198 

2,726 

13 

3,099 

0,803 

8,447 

2,726 

14 

3,594 

0,803 

9,797 

2,726 

15 

4,126 

0,803 

11,246 

2,726 


From the illustration of the three waves, for a wave with 
large wave amplitude, H = 2A can not be confirmed 

although j = ’ lmax Vmin < 2 could happen and although 
the profile is asymmetrical. 


III. WATER WAVE-SURFACE EQUATION 
In this section, the formulation of water wave surface is a 
revision of the formulation of Hutahaean (2019a), where 
some typographical errors occurred. The water wave 
surface equation is obtained by integrating KFSBC, 

3 r] Gk 

— =- fi^^yqjacoskxsinat 

— asinkxsinat —.( 1) 

yo dx 

/?!&)) = sinhkth + tj).(2) 

fi(j]) = coshk Qi + rf) .(3) 

yis weighting coefficient at weighted total acceleration 
equation with a value around 2.202-3.0. This research uses 
y = 2.483 (Hutahaean (2019c,d)). 

(1) was integrated against time-t where acccording to wave 
number conservation, i.e. at z =p, = 0(Huatahean 

(2019a)), which means that /? 1 (t))and^(t))have constant 
values against time t, thus the integration of the first tenn 
right side of the equation was completed by integrating 
sinat. 


r)(x,t ) = — fi 1 (j]')coskxcosat 
ya 

The integration of the second termright side of the equation 
will be completed using partial integration method, as 
follows 


Assume a function of f = cosat—. This cosertfunction 
1 ax 

was used so when it was differentiated against time t, 
sinat — will be formed, i.e. 

3x 

a / a?7 a 2 ?? 

— = —asinat -h cosat - 

3f 3x 3t3x 

This differential equation was multiplied with dtand 
integrated against time t, 

r 377 r 3 2? 7 

f = —a sinat — dt + cosat - dt 

J 3 x J 3t3x 

Substitute /and the first term right side was moved to the 

left and /was moved to the right and both equations were 

divided by a, 

f 3r/ 1 3p If 3 2 rj 

sinat — dt8 =- cosat -1— cosat - dt 

J 3x a 3 x a) 3t3x 

The integration of the second termright side of the equation 

can be completed the same way, but with an assumption 


that 


at 2 ape 


is a very small number, the integration can be 


completed by integrating just th ecosat element. 

3rj 1 3p 1 3 2 p 


/ 


sinat —-dt 
3x 


■ cosat —- H —-sinat - 
a 3x a 2 3t3x 


Substitute the result of integration, 

, , Gk 

rjKx.t) = —coskxcosat 
ya 


Gk 

ya 1 


-fi(j]')sinkx 


3p 1 3 2 r] \ 

-cosat -1— sinat - 

3x a 3t3x J 


Working on an assumption that -sinat — 4s a very small 

o 3t3% 

number and can be ignored, 

Gfc Gk 3rj 

rjbc,t) = —fii \jj) coskxcosat H- fiujjsinkxcosat — 

The equation was differentiated against horizontal-x axis. 

3r) Gk 

— =- fiArjjksinkxcosat 

3x ya 

Gk 3p 

H- Bynjkcoskxcosat — 

ya 3x 

Gk 3 2 r] 

H- B (ji) sinkxcosat —- 

1 3 x 2 


or 

377 

3x 


— — fi^rj^ksinkx + — fi(.rj)kcoskx^- + 
— /?&))sinkx jj-j) cosat .(4) 


ya 


In accordance with the provision at the velocity potential 
equation where there is tfunction only,* function 
onlyandzfunction only, then at the water wave surface using 
variable from velocity potential equation, water surface 
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equation also has variable that is tfunction only anxfunction 
only, so it can be stated that the general forniof (4) is, 
a?7 

— = fix)cosat 
3x 

In this equation fix) is just a function of x. fix) in (4) is, 
fix) = 

(Gk Gk 

- /3 1 ir])ksinkx H- pig)kcoskxcosatf ix) 

V ya ya 

ya 3x> 

Bearing in mind p 1 irj) and pig )should be constant against 
time land against horizontal x axis. In the second and third 
terms right side of the equation there is a function of time t, 
i .e.cosat , then the term shouldn’t be there. 
Therefore/Oc) is, 

Gk 

fix) -- [f irfksinkx 


ya 


Thereby 
3g Gk 


Pi irf) ksinkxcosat 


3x ya 
Substitute this equation to (1), 
3 rj Gk 


at 


ya 


Pi irf) acoskxsinat 


+ 


(Gky 

(—) pig)Piig)aksin 2 kxsinatcosat 
\yo 1 


is integrated against time t. Integration will be completed 
by bearing in mind wave number conservation , thus 
integration was done by integrating only sinusoidal element 

and bearing in mind that sinatcosat = ^ sirilat , 

Gk 

n ix, t) = —/?, in) coskxcosat 
ya 

1 /Gfc\ 2 


4 \ya / 


P irf)Pi irf) ksi n 2 kxcos2at 


..(5) 


To determine the values of Piirf), pirf), an approach was 
done that water surface equation is sinusoidal, i.e. 
r) 0 ix,t) = Acoskxcosat .(6) 

Therefore, the approach of hyperbolic function is Piirf) = 

Piiri 0 ), Piv) = PiiJo)- 

Gk 

T)ix, t) — — Piir) 0 )coskxcosat 


ya 

2 


) Piri 0 )Piiri 0 )ksin 


2 kxcos2at 


..(7) 


IV. THE CALCULATION OL G, k AND A AT 
BREAKER POINT 

The calculation of G, k andAat breaker point was done 
using shoaling-breaking model (Hutahaean (2019b)). In 
principle, the model is a transformation analysis of G, 
fcandAfrom the deep water to breaker point. 

4.1. The calculation of Gandfein Deep Water 
To calculate G and fctwo equations were needed. KFSBC 
and surface momentum equation, with the formulation can 
be seen in Hutahaean (2019b). 


Wave number kin the deep water was calculated using the 
following equation, 

V a = gk .(8) 

gis gravitation velocity, whereas Ais wave amplitude. This 
equation is a quadratic equation of wave number 
kwherekcan be calculated with a simple equation, i.e. 
finding the root of the equation. Wave amplitudemaximum 
at (8) is (Hutahaean (2019a)) 


0-9 Iff 
2y 2 er 2 


( 9 ) 


Therefore, wave amplitude in (8) should be smaller or less 
than A max . 


Furthermore, deep water depth should be determined which 
is the starting point of shoaling-breaking analysis where 
wave amplitude Ais known as input of the model. Deep 
water depth is calculated with, 
h 0 = h y L = 1-65 nA max .(10) 

The coefficientof 1.65 was obtained from the calibration of 
the breaker depth of the shoaling-breaking model against 
breaker depth of SPM (1984), Hutahaean (2019b). 

After wave number k was obtained, wave constantGwas 
calculated in the deep water with the following equation 

yaGp g) = gA .(11) 

Where wave amplitude A is known. 


4.2. Shoaling-Breaking Model Equations 
In the shoaling-breaking analysis, there are 3 (three) 
variables that change along with the change in the depth, 
i.e.k, AandGwhere there is a dependency among the 
changes of the three variables. 

To calculate the three variables, three conservation 
equations were used, and the formulation of shoaling- 
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breaking equations in this section can be seen in Hutahaean 
(2019b). 


The first equation is the result of the substitution of KFSBC 
and energy conservation equation to wave number 


conservation equation, i.e.: 


( 

^ 7l\ Gkf\dk 3 h 

(h + - +- —= -k — + 

\iGk 3 f 

\ 

.(12) 

2/ 2per J 5x 5x 

2ya 

.(13) 

r-(*&)- « + 

kA\ 

2 ) 


(12) and (13) were formulated using velocity potential at 
the sloping bottom, where in this case, 



—is bottom slope with negative value. After—is knownthe 

3% 3x 


second equation is done, i.e. energy conservation equation 

£=-» Gk .^ 


b = 



■( 18 ) 


After—-and—were known — can be calculated using the 
ax ax ax b 

third equation, i.e. KFSBC, 


3A /3G , _ 3 k\ r 

y°—= (—fc + c—)/ 

3x V3x 3x / 


.(19) 


The calculation of—, with(12), was done iteratively .At the 
beginning of the calculation, pin (18), is not known since 
^-has not been known, as initial approximation 
equation—= — -—was used andpwas calculated with (18) 
and by obtaining the new —from (12), pwas recalculated 

3x 

with (18). This calculation is repeated until a stable —is 

obtained where, in general, stability is achieved after 5 

times iteration. After stable— anduare obtained, —is 

ax b ax 

3/4 

calculated with (17) and—with (19). Furthermore, wave 

3x 


number, wave constant and wave amplitude are calculated 
at the point x + &cusing Taylor series, i.e. 

a k 

k x+Sx = k x + 8x — 

5G 

G x+Sx = G x + 8x — 

a A 

A x + Sx = A x + Sx — 

V. THE RESULT OF THE MODEL 

5.1 Illustration of the Result of Shoaling-Breaking Model 
As an illustration of the result of shoaling-breaking model, 
shoaling and breaking analysis was done for the wave with 
wave period of 8 sec, with deep water wave amplitude A 0 = 

1.056 m, obtained from (6). Sea bed slope — = —0.005,the 

dx 

result of the model is presented in Fig.3., where at the 
beginning there is shoaling, i.e. the enlargement of wave - 
amplitude A, followed by breaking. Output of the model at 
breaker point includes: breaker amplitude A b = 1.402 m, 
breaker depth h b =3.15 m, breaker wave number k b = 
1.4355 , wave constant G b = 0.00139. 


1.5 

E. 

< i 

aT 

■a 

0.5 

a. 

E 

< A 

















QJ ^ 

| 0 5 10 15 

Water depth, h (m) 


Fig. 3. Wave amplitude as the result of shoaling-breaking 

model 

5.2. The Calculation of Breaker Height H b 
As has been stated before, the output of shoaling-breaking 
model is breaker wave amplitude A h that must be changed 
into breaker height H b . 

The calculation of H b was done using water wave surface 

equation (1). With wave parameter input at breaker point, 

i.e. A b , h b , k b andG b , the elevation of wave crest t) max , the 

elevation of wave trough r] min were calculated, then// h = 

Vmax ~ r lmin- From the example of the result of shoaling- 

breaking in fig.3, breaker height H b = 2.522 m was 

obtained, with the value of—= 1.798, where//,, =£ 

Ab 

2A b . The wave profile at breaker point is presented in Fig.4. 
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which shows the almost perfect cnoidal shaped of wave 
profile. 



Fig. 4. Wave profile at breaker point 

The value of^ = 1.798is quite close to H b = 2A b . 

Therefore in Hutahaean (2019b), H b model was found to be 
quite close with the breaker from empirical equation 
formulated from laboratory observation. 

5.2. Calibration and Adjustment of the Result of the Model 
As a comparison of the result of the model, breaker height 
index equations from 5 (five) researchers were used and for 
the comparator breaker height, the average value of breaker 
height produced by the five breaker height index was used. 
The breaker height index (BHI) equations that were used 
are equations from Komar and Gaughan (1972), Larson, M. 
and Kraus, N.C. (1989), Smith and Kraus (1990), Gourlay 
(1992) and Rattana Pitikon and Shibayama (2000). As a 
comparator of the result of breaker depth, equation from 
SPM (1984) was used. 

The adjustment of breaker height from the result of the 

model with the result of breaker height index equation was 

done by changing the value of wave constantG h , and that in 

order for the breaker height to be correspond to breaker 

height from breaker height index, G b from the result of 

shoaling-breaking model must be multiplied by 0.76. 

Whereas breaker depth of the output of the shoaling- 

breaking model was also determined by the choice of deep 

water depth or the initial depth that was used. The result of 

the calibration of breaker depth model with breaker height 

fomi SPM (1984) shows that by using deep water depth 

h 0 = ——breaker depth that corresponds to SPM (1984) 
k 0 

was obtained. 

5.3. Some of the Model Results 


In this section, model was done for several wave periods, 
with deep water wave amplitude from (6), and the result of 
the calculation is in Table (5) and Table (6). Adjustment 
with the result of breaker height index (BHI) equation was 
done with a process that has been discussed in sub-chapter 
(5.2), i.e. by multiplying wave constant G b with 0.76. The 
result is that breaker height model is very close with breaker 
height from BHI, where breaker height from BHI is the 
average values of 5 breaker height index equations as stated 
in (5.2). 

Whereas breaker depth, the result of the model is very close 
to breaker depth from SPM (1984), that was obtained using 
deep water depth h 0 = 

k 0 

Table (6) shows that the model produces the value — = 
0.801for all wave period, which is very close to the result 
of BHI equation, i.e. ~ L= 0.802. 


Table.5: The Result of Shoaling-breaking Calculation 


T 

^0 


(m) 

h b 

(m) 

(sec.) 

(m) 

model 

BHI 

Model 

SPM 

7 

0,809 

1,93 

1,929 

2,411 

2,406 

8 

1,056 

2,522 

2,52 

3,15 

3,142 

9 

1,337 

3,194 

3,19 

3,988 

3,977 

10 

1,65 

3,943 

3,938 

4,922 

4,909 

11 

1,997 

4,77 

4,765 

5,956 

5,94 

12 

2,377 

5,675 

5,67 

7,09 

7,07 

13 

2,789 

6,659 

6,655 

8,32 

8,297 

14 

3,235 

7,722 

7,718 

9,65 

9,622 

15 

3,714 

8,862 

8,86 

11,08 

11,046 


Table.6: The Result of Shoaling-breaking Calculation 


(continued) 


T 

(sec.) 

Ab 

(m) 

Hb 

K 


H b 

Model 

BHI 

L b 

Ab 

7 

1,073 

0,801 

0,802 

0,569 

1,798 

8 

1,402 

0,801 

0,802 

0,569 

1,799 

9 

1,775 

0,801 

0,802 

0,569 

1,8 

10 

2,191 

0,801 

0,802 

0,569 

1,8 

11 

2,651 

0,801 

0,802 

0,569 

1,799 

12 

3,155 

0,8 

0,802 

0,569 

1,799 

13 

3,702 

0,8 

0,802 

0,569 

1,799 

14 

4,294 

0,8 

0,802 

0,569 

1,798 

15 

4,929 

0,8 

0,802 

0,568 

1,798 
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ttt 11 . Hh _ , „ 1.788 

Wave steepness at breaker point is—= 0.569 =-, 

which is very close to or similar to the analytical result as 
shown as follows. In (10) there is a breaking characteristic 
when: 



/?! (y^)and P (yy)i s t ^ le va l ue i n the deep water 


needed in shoaling-breaking model as well as water wave 
surface equation. 

The calculation of wave force in a stmcture, should take 
into account the cnoidal-shaped wave profile, so that a 
better estimation on the location point of wave force is 
needed. Furthermore, the planning of coastal constmction 
elevation should take into account the asymmetry wave 
profile so that a parameter as big as half the wave height 
cannot be used. 


where- 




tanhk 


(h + y) = 1. then breaking occurs 


when 


kA 1 

T “ 1 + p 


or 

A„ _ 1 

L b (l + p) 7T 

The value — = 1.80, obtained 

Ab 

H„ 1.80 

L b (l + p) 7T 

With an assumption that the bottom slope is very small, 

H„ 1.80 
— «-= 0.572 

L b IT 

It shows that the numerical result is very close to or similar 
to the analytical result which also proved thatj^- = 1.80is a 


good value. 


VI VI CONCLUSION 

Wave profile at breaker point is cnoidal-typed wave, where 
wave crest is not symmetrical with wave trough. 
Interpreting the result of shoaling-breaking model should 
perform analysis on water wave surface profile. As an 
approximation^ criteriathat at breaker point the 
relation—= 1.8 applies with wave steepness — = 0.572 

A b L b 

can be used. 

Water wave surface equation that was used can produce 
wave types that corresponds to the Wilson criteria (1963), 
and compatible enough with the shoaling-breaking model 
that was used. However, further research is still needed 
because adjustment should still be done on wave contant G 
of the shoaling breaking model. Further research is still 
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